Measurement of pd — > 3 He rj in the Sn Resonance 
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Abstract 

We have measured the reaction pd — > 3 Her/ at a proton beam energy of 980 
MeV, which is 88.5 MeV above threshold using the new "germanium wall" detector 
system. A missing-mass resolution of the detector system of 2.6% was achieved. 
The angular distribution of the meson is forward peaked. We found a total cross 
section of (573 ± 83 (stat.) ± 69 (syst.) ) nb. The excitation function for the present 
reaction is described by a Breit Wigner form with parameters from photoproduction. 
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The production of 77-mesons is interesting because it opens the possibility of studying 
the interaction between the lightest isoscalar particle and the nuclear environment. Haider 
and Liu were the first to show that even bound 77-nucleus systems, i.e. ?7-mesic nuclei, 
could be possible [|TJ . Based on the results of Bhalerao and Liu || they found an attractive 
77-N interaction which in their calculations leads to bound states for nuclei with mass 
number A > 10 [[|. Rakityanski et al. [Q even relaxed this condition to A > 2. The 
widths of such states were predicted to be narrow enough to be observable for nuclei with 
A > 4. Wycech et al. || also predicted the formation of mesic nuclei in dd — > 4 He rj, but 
not in pd — » 3 He?7. In contrast, Abaev and Nefkens ||, as well as Wilkin |7J showed that 
the formation of quasi-bound ??- 3 He states in the reaction pd — > 3 He r] should indeed be 
possible. 

In addition, the reaction pd — > 3 He 77 is of interest due to its surprisingly large cross 
section close to threshold making this reaction a prime candidate for the source of 77- 
mesons in tagged 77 facilities ||. 

A detector system called the "germanium wall" |J was built at the COSY facility in 
Julich (see Figure ). In its complete setup, the germanium wall is a stack of four position 




Figure 1: The Germanium Wall The detector system "germanium wall". In the present 
measurement detector E2 was removed. 



sensitive high-purity germanium detectors having a conical acceptance with an opening 
angle of ±287.5 mrad. In the centre of each detector is a hole with a size of ±28 mrad 
allowing the primary beam to pass through. Two types of detectors are used, one 1.3 
mm thin diode ( "quirl-detector" ) for determining the reaction vertices through its good 
position resolution given by the crossing of two counterrotating spirals and three 17 mm 
thick diodes for measuring the particle energies ("energy-detectors"). For further details 
see Ref . f| . The setup used for the present measurement consisted of one quirl and two 
energy-detectors (Quirl, El and E3, see Figure ). 

First measurements with the "germanium wall" showed the good missing-mass res- 
olution of the system. The reaction pd — > 3 He?7 was studied at a proton beam energy 
of 980 MeV (88.5 MeV above threshold) leading to almost Att acceptance of the detector 
system for the product 3 He-particles. We performed two runs at different times. 

The target was a cell filled with liquid deuterium with 6 mm diameter and thicknesses 
of 2.4±0.2 mm (run A) and 4.4±0.2 mm (run B), respectively |TIJ. The COSY extracted 
proton beam was focussed onto the target yielding a spot with a radius a = 0.5 mm 
and a divergence of 6 mrad. These parameters together with the short distance between 
target to detector yields a total angular uncertainty of 16 mrad, where the individual 
contributions are linearly added. This uncertainty is much larger than that resulting 
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Figure 2: Particle identification using a AE-E spectrum measured with the first two detectors 
of the "germanium wall". Protons, deuterons, tritons and 3 He-particles can be identified as 
indicated in the figure. 



from the position resolution of the detector, which is in the order of 2 mrad. The beam 
had a momentum spread of Ap/p = 8 x 10 -4 [TT| . 

The energy and direction of the emerging 3 He-particles were measured by the "ger- 
manium wall" . Figure [2] shows a AE-E spectrum demonstrating the capability of the 
detector system for particle identification. Through energy and emission direction mea- 
surement of 3 He-particles, the missing-mass was calculated. A missing-mass spectrum 
for run B is shown in Figure |3|. The 77-peak is clearly visible with a resolution of 
a = (6.1 ± 0.5)MeV/c 2 . Background is mainly caused by multi-pion production (e.g. 
pd — > 3 He7r + 7r~, pd — > 3 He7r°7r°, etc.). Low-energy 3 He background events at small 
angles were not detected, because of the minimum opening of the detector system. These 
events correspond to large relative energies between the pions and thus to large missing- 
mass values. Therefore, the spectrum is truncated at 600 MeV/c 2 . For run A a missing- 
mass resolution of a = (8.2 ± 0.6)MeV/c 2 is obtained. The whole body of data was 
divided into 5 and 6 angular bins for run A and run B, respectively. For each bin a Gaus- 
sian together with a background function were fitted to the corresponding missing-mass 
spectrum above 450 MeV/c 2 . For run B, two different shapes for the background were as- 



sumed: a polynomial of third order and a function BG 
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mm the missing-mass and <2j parameters to be fitted. Both functions lead to the same 
results. For run A, the background is not so clearly separated for large missing-masses 
from the rj peak as it is in run B. Therefore, several functions were tested. Polynomials 
were fitted to the range lower than the i] peak and only the first point above the peak. 
Alternatively a step like function (cummulative Lorentzian) was fitted to all data. The 
r) peak was always assumed to be a Gaussian. Finally the number of pd — > 3 Her] events 
was obtained by integrating the Gaussians and weighted means were deduced. Further 
details of the data analysis procedure are given elsewhere \T2\. 
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Figure 3: Missing-mass spectrum. The 77-peak at approx. 547 MeV/c 2 is visible with a 
resolution of (14.2 ± 1.2) MeV/c 2 FWHM. The background is caused by multi-pion production. 



Due to beam halo during the experiments the intensity of the beam had to be reduced 
to a level of 10 5 protons per second. Thus pile up and detector damage were avoided, 
but event statistics was strongly reduced. The measured angular distribution is shown in 
Figure [|. The error bars shown represent the statistical errors only. In addition, there are 
systematic uncertainties: target thickness 10% and 5% in the different runs, respectively, 
luminosity calibration 7%, corrections due to trigger and detector inefficiencies due to 
nuclear interactions in Germanium (see Ref. ||13|1 ) 5%. The total systematic error of 13% 
and 10% (when added in quadrature) for the two runs, respectively is smaller than the 
statistical error. The efficiency of the data analysis (80%) was studied by Monte Carlo 
simulations. The simulated detector response was found to be in excellent agreement with 
the experiment Jl4|. Before and after each run the detector was checked with radioactive 



sources. No significant deviations in the amplifier gains were found as expected, since the 
electronic circuits were kept at a constant temperature. Since both runs were performed 
under different experimental conditions with respect to beam halo, target thickness and 
distance between target and detector |L2| , different systematic errors lead to an enhance- 



ment of run A compared to run B which is slightly above the statistical error given above. 
The halo was 2.2 times more intense during run A than in run B thus leading to a larger 
combinatorial background than in run B and hence to larger error bars of the integrated 
Gaussians for the rj peak. 

The angular distribution is forward peaked. This is in contrast to the results of Mayer 
et al. m reporting almost isotropic distributions close to threshold. The cross sections of 
Banaigs et al. |15| measured through dp — > 3 He?7 at a slightly higher excitation energy 
(corresponding to an equivalent proton kinetic energy of 1047 MeV) agree with our data. 
When we fit Legendre polynomials to the present data, an unphysical negative value 
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Figure 4: Angular distribution for the reaction pd — » 3 He i] at 88.5 MeV above threshold in 
the centre of mass system. The solid curve shows a Legendre polynomial fit to the data, the 
dashed curve a calculation within a resonance model (see text). 



for cos(6 l ) = — 1 is obtained. In order to overcome this deficiency we have added one 



data point at cos (#)=-! from Ref. ||T6[, who had measured an excitation function for 



this angle. This point is also shown in Fig. |j. In order to take into account the fact 
that the point was measured at a slightly different energy of 0.5%, we have doubled its 
statistical error. Another data point in the literature taken at a somewhat higher beam 
energy (see Ref. |Io|) is also shown in Fig. [|. Again its error bar was doubled. A Legendre 
polynomial fit to all points of the angular distribution, i.e. the present ones from both runs 
as well as the one from Ref.'s [|l^, [L6|] weighted by their total errors, yielded parameters 
A = 45.6 ± 5.9, Ai = 47.3 ± 10.9, and A 2 = 5.8 ± 8.5 all given in nb/sr, and with a 
X 2 /nfree = 0.3 from which a total cross section of (573 ± 74) nb follows. In addition to 
this statistical error the systematic error is assumed to be 69 nb. The result is insensitive 
to the added point, because the differential cross section is small for the backward angle 
emission. Including higher degrees in the fit procedure does not improve the fit, a lower 
degree gives an unphysical negative value for cos(6>) = — 1. 

Kingler |^0| has extended a model originally developed for pion production |2l| to 
include higher-resonances. The original model was limited to only pion exchange, A 
resonance excitation as well as non resonant contributions. The extension treats other 
meson exchanges as well as higher nucleon resonances than the A. The vertex function 
for the different baryon-baryon-meson couplings were calculated in a simple quark model 
and are momentum dependent. For the present reaction the largest contribution to the 
cross section comes from the NNp and NNu interactions while the contribution due to 
NN*(1535)ir interaction is one order of magnitude smaller. The contributions of other 
resonances like A^*(1440), A r *(1650), and A r *(1710) are even smaller. However, the form 
factor is calculated only for harmonic oscillator wave functions with the frequency being 
a free parameter varied to fit the experimental data. The model predictions are shown in 
Fig. |] as dashed curve. Obviously, the calculation shows structures not observed in the 
present data. 

The extracted total cross section together with earlier data from Mayer et al. || is 
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shown in Figure [5[ From the angular distributions given by Banaigs et al. ||15|| , Loireleux 
T?j and Kirchner |18|] we extracted total cross sections by fitting Legendre polynomials. 



The results are also shown in Figure [|. Kirchner claimed that the data from Ref. flT7| 
suffer from important electronic problems. No further details are given. The data from 
Mayer et al. || indicate a strongly rising cross section close to threshold. 

Also shown in the Figure is a normalized calculation within a two-step model de- 
veloped by Kilian and Nann [IS]. Within this model a pion is produced in a first step 
through pp — > dir + . In a second step, this pion produces the r/ in an interaction with the 
neutron. A kinematical velocity matching yields the maximum close to threshold. The 
present data do not support this model as the dominant reaction mechanism. 




Figure 5: Total cross sections for pd — ► 3 Her/ reaction as function of the proton beam energy. 
Data are indicated by different symbols. The kinematical velocity matching in a two step process 
is shown as dotted curve. The solid curve is a calculation employing the matrix element from 
photoproduction on the proton. 

The energy region from 900-1100 MeV corresponds to the centre of the N* Sn res- 
onance (r ~ 200 MeV) known to couple strongly to the r\- N channel ||22|| . One may 
therefore attempt to describe the cross section by an intermediate iV*(1535) resonance 
excitation. The cross section is calculated as 

a(E) = ^\M(E)\ 2 (1) 
P P 

with E the excitation energy and M the matrix element. All momenta p are in the centre 
of mass system. This is calculated as in photoproduction on the proton [ p3[ 

AT 2 

|M(£)|2 = (E- m ,Av(Ef V 

with 

Y{E) = T R (b v ^ + + bj\ . (3) 

V Pr,,R P*,R J 

Similar to Ref. |23| , we applied a width at the resonance of T#=200 MeV, a Breit- Wigner 
mass of 771^=1540 MeV/c 2 . The branching ratios were set to 6^=0.47 for the rj decay, 
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6,r=0.48 for the pion decay and 6 ff7T =0.05 for the two pion decay [521. The momenta at the 
resonance position are indicated by the index R. The only free parameter is the strength 
A taken to be 241 nb in order to fit to the present data point. The calculation is shown 
in Figure |5] as a solid curve. The trend of the data is reproduced, which may be taken 
as an indication that production of the iV*(1535) resonance is the dominant reaction 
mechanism and that the product of kinematics and form factor changes only very little 
over the present energy range. The overall shape of the calculation slightly underestimates 
the data of Mayer et al. || close to threshold. For a more detailed investigation of fine 
structure, additional data in this region is needed. An enhancement close to threshold 



was also seen in r] production in NN interactions |24|,|25| and was attributed to a strong 
final state interaction. The agreement between the excitation functions for pd — > 3 B.erj 
and 7p — > pr] reactions excludes strong FSI between the nucleus and the rj except for the 
near threshold region. 
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